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Abstract 

> ; 

The understanding of the large-scale structure formation requires the resolution of coupled 
nonlinear equations describing the cosmic density and velocity fields. This is a complicated problem 
that, for the last decade, has been essentially addressed with N-body simulations. There is however 
a regime, the so-called quasi-linear regime, for which the relative density fluctuations are on average 
^sO , below unity. It is then possible to apply Perturbation Theory techniques where the perturbation 

expansions are made with respect to the initial fluctuations. I review here the major results that 
have been obtained in this regime. 

O ' 1 The Gravitational Instability Scenarios 

a: 

Perturbation Theory (PT) has had important developments in the last few years simultaneously from 
theoretical, numerical and observational points of view. These techniques allow indeed to explore 
nonlinear features revealed in the statistical properties of the large-scale structures of the Universe, 
for which precise and robust data are now available. The general frame of these calculations is based 
on the dynamics of a self-gravitating pressure-less fluid. The large-scale structures are then assumed to 
have gravitationally grown from small initial fluctuations. It is important to note that in the following 
these initial fluctuations will be assumed to follow a Gaussian statistics. That excludes a priori exotic 
models that make intervene topological defects as seeds of structure formation. 

By linearizing the field equations, i.e. when the quadratic coupling between the fields is neglected, 
one can compute the growth rate of the fluctuations. This allows to compare for instance the Cosmic 
Microwave Background temperature anisotropies as measured by the COBE satellite with the local 
density fluctuations observed in galaxy catalogues. The linear approximation cannot be used however 
to explore the statistical properties of the density field: the local fluctuations are just amplified, their 
shapes are not changed, and therefore the linear density field remains Gaussian if it obeyed such 
statistics initially. 

2 The Perturbation Theory 



The principles of these calculations have been initially presented by Peebles [35], then explored in 



more details in a series of recent papers |22|, [L2], |||, [27|, |28|, ^3|, ||, [37|, 35]. The starting point of all 
these calculations is the system of field equations, Continuity, Euler and Poisson equations, describing 
a single stream pressure-less fluid. 



The density and velocity fields which satisfy those equations are then expanded with respect to 
the initial fluctuation field, 



*(t,x) = £>«(i,x), u(t,x) =$>W(t,x). (1) 

The fields S^(t, x) and uW(t,x) are just the local linearized density and velocity fields. They are 
linear in the initial density field. The higher order terms, 5^ 2 \ 5^ 3 \.. are respectively quadratic, cubic,., 
in the initial density field, and, therefore, do not obey a Gaussian statistics. 

2.1 The Density Field 

The time and space dependences of the linearized density field, factorize so that it can be written, 

tfW(i,x) = D(t) J d 3 k5(k) exp(ikx), (2) 

where 5(k) are the Fourier transforms of the initial density field. They are assumed to form a set 
of Gaussian variables. Their statistical properties are then entirely determined by the shape of the 
power spectrum, P(k), defined by, 

(5(k)5(k')) =5 Dhac (k + k')P(k), (3) 

where (.) denotes ensemble averages over the initial conditions. The function D{t) is determined by 
the cosmological parameters. It is proportional to the expansion factor a{t) for an Einstein-de Sitter 
universe. In general it has been found to depend on the cosmological density Vt in such a way that 
(dlog£>/dloga) ~ 6 H and to be very weakly dependent on the cosmological constant A |2J| . 

The higher order terms can all be recursively obtained from the linear term. In general one can 
write the i th order term as [E2|, 



*W(t,x) = £>*(*) J d 3 k 1 ...d 3 k i <5(k 1 )...<5(k i ) exp[ix-(k 1 + ... + k i )]F(k 1 ,...,k i ), (4) 

where F(ki, . . . , kj) is a homogeneous function of the angles between the different wave vectors. Note 
that the time dependence given here is only approximate for a non Einstein-de Sitter Universe. It is 
however a very good approximation [||, |l2], |4|, [B| . 

2.2 The Velocity Field 

For the velocity field the situation is very similar. We have first to notice that in the single flow 
approximation the vorticity is expected to be diluted by the expansion (eg. |35|]) and thus to be 
negligible at any order of the perturbation expansion. Then it is more natural to present the properties 
of the velocity field in terms of the local divergence (expressed in units of the Hubble constant), 

_ V x • u(t,x) 

^' X)= H(t) ■ (5) 



We then have, 



d log a 



The higher order terms can be written 

#«(i,x)«tt°- 6 (i) D\t) J d^.-.d^kx)...^) exp[ix-(k 1 + ... + k,)]G(k 1 ,...,k i ), (7) 

where G is another homogeneous function, different from F. Note that in general the time dependence 
of 0® is not the power i of the time dependence of 9^ . We will see that it induces remarkable statistical 
properties for the local velocity field. 



2.3 Implications for the Statistical Properties of the Cosmic Fields 



In general the consequences of the existence of higher order terms can be separated in two categories: 
these terms affect the mean growth rate of the fluctuations and introduce new statistical properties 
because of their non-Gaussian nature. I will briefly review both aspects here. 

3 The Mean Growth Rate of the Fluctuations 

The mean growth rate of the fluctuation can be calculated from the shape and magnitude of the 



variance of the local density contrast [26, 32, 38], or equivalently of the evolved power spectrum 



1 33, 25]. For instance the variance, (5 2 ), can be calculated using the expansion in (p]), 

(tf 2 ) = <(*(*) +*( a > + ...) 2 >. (8) 
By re-expanding this expression with respect to the initial density field one gets, 

( 5 2 ) = ( ( 5 (i)) 2 )+ 2 ($« ^) + ((V 2 )) 2 >+ 2{(5W^)V- (9) 

The term (5^ 6^') is zero because it is cubic in the initial Gaussian variables. The calculation of the 

corrective terms (the so-called "loop corrections", see [22, 37, ||]), ( (<^) ) + 2(<5 ( - 1 * ) 5^}, can be 
done from the recursive solution of the field equations. In practice these calculations become rapidly 
very complicated but they can be simplified in case of a power law spectrum, 

P(fc) oc k n . (10) 

Thus for n = —2 one gets 

{d2) = {dL) + iH ( ^ )2 w ( ^ n> + °- 88 {6L)2 (n) 

where 5n n is the linearized density contrast, <5n n = This prediction can be compared with nu- 

merical results obtained from iV-body codes. An interesting way to compared the two is to use the 
phenomenological description of the nonlinear growth rate of the fluctuations proposed by Hamilton 
et al. (2^]. This description relates the mean linear density contrast at the scale R\ m to the full 
non-linear one at the scale i? non _ii n defined by 

^non-lin (l + {?)) = #lhr (12) 

The universal transform proposed by Hamilton et al. (given here in terms of the power spectrum) is 
presented in Fig. 1 (solid line) and compared to the prediction (|TT]), (short dashed line). The latter 
gives remarkably well the position of the departure from the pure linear regime (straight line). On the 
other hand the position of this transition is quite poorly given by the Zel'dovich approximation (long 
dashed line). 
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Figure 1. Comparison of the PT predictions (|ll]) with the Hamilton et al. p^] prescription for the 

growth rate of the fluctuation (solid line) in case of P(k) oc k~ 2 . The short dashed line is the 
prediction of the next-to- leading order Perturbation Theory result |38| , and the straight line is the 

linear theory prediction (figure taken from |38|]). 

Surprisingly the corrective terms are finite only for n < — 1. For n larger than —1 the loop terms 
contain a divergence so that the result explicitly depends on a cutoff k c introduced in the shape of the 
power spectrum at large k. Thus when (R k c ) is large we have, 

(S 2 ) = (tL)+C n (k c R) n+1 {5lf (13) 

where R is the filtering radius. The existence of this divergence was not really expected since the 
numerical results do not indicate any significant change of behavior for n ~ — 1. One could argue that 
in practice this is not a relevant problem since in realistic scenarios (like the CDM power spectrum) 
the index reaches n = —3 at small scale, so that the corrective terms are naturally regularized. 
However with such an interpretation it implies that the growth of the very large scale fluctuations is 
fed with smaller scale fluctuations, which would be slightly surprising in view of the numerical results 
accumulated over the last few years. Another possible explanation is that the higher order corrections 
induce other divergences that cancel each other (a phenomenon quite common in statistical physics). 
Then the corrective terms to the linear growth rate would not be proportional to the square of (<5f in ) 
but to a smaller power of it. Future analytic investigations may be able to throw light on this problem. 

4 The Emergence of non- Gaussian Features 

4.1 The Moments 

The other major consequence of the existence of non-linear terms in (||) is the apparition of non- 
Gaussian properties. Although it is possible to characterize non-Gaussian features in many different 
ways, most of the efforts have been devoted to properties of the one-point probability distribution 



function (PDF) of the local density. More particularly, attention has been focussed on the moments 
if this distribution and how they are sensitive to non-linear corrections. In his treatise [35|, Peebles 
already considered the implications of second order perturbation theory for the behavior of the third 
moment of the local density in case of initial Gaussian fluctuations. Indeed, using the expansion ([]]), 
the third moment reads, 

( 5 3) = ( (jd) + 6 (2) + . . .) 3 ) = ( (^)) 3 ) + 3 ( (^)) 2 5^) + ..., (14) 

and what should be the dominant term of this series, (CsQ-A ), is identically zero in case of Gaussian 
initial conditions. As a consequence the third moment is actually given by 

(5 3 ) « 3((>))V)}. (15) 

Peebles calculated this expression neglecting the effects of smoothing and found (for an Einstein-de 
Sitter Universe), 

<* 3 > - y (C) 2 , (16) 

so that the ratio, (5 3 }/(5 2 ) 2 , is expected to be finite at large scale. It is actually possible to extend 
this quantitative behavior to higher order moment and to show that the cumulants (i.e., the connected 
parts of the moments) are related to the second moment so that the ratios, 

S P = (S p )J{S 2 ) p -\ (17) 

are all finite at large scale. As mentioned before, the coefficient 53 was computed by Peebles p5|] , Fry 
flq| derived S4 and eventually Bernardeau || gave the whole series of these coefficients. 

Unfortunately these early calculations did not take into account the filtering effects, that is that 
the ensemble averages should be done on the local smoothed fields. This problem was addressed 
numerically by Goroff et al. [^] for a Gaussian window function for the third and fourth moments. 
More recently these two coefficients have been calculated analytically and semi-analytically in [28, 31] 
for this window function. However, the results turn out to be simpler in case of a top-hat window, 



as it was noticed by Juszkiewicz et al. [27| for S3 and power law spectra. The coefficients S3 and 
S4 were calculated for this window in Q for any spectrum and any cosmological models. Eventually 
Bernardeau || proposed a method to derive the whole series of these coefficients from the spherical 
collapse dynamics. 

I recall here the expression of the first two coefficients S3 and S4, as a function of the shape of the 
second moment, 

s 3 (R) = y + (18) 

60712 62 7 (if) 7 7 (R) 2 2 d 7 (R) 
Si{R) = 1323" + + + 3 dlogl?' (19) 

with 

One can notice that S3 depends only on the local slope, and that S4 also depends (but weakly) on the 
variations of that slope. These results are valid for an Einstein-de Sitter Universe but are found to be 
weakly Q and A dependent |L2], §, |] • 
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Figure 2. Comparison of PT predictions with results of numerical simulations with CDM initial 
conditions for moments of the local density contrast. The PT results (dashed lines) are compared to 
the measured S p coefficients as a function of radius for three different time-steps (figure taken from 

§)• 

These coefficients have been tested against numerical results by different groups ||, ^, |3l| and 
have been found to be in extremely good agreement with the numerical results. In particular as shown 
in Fig. 2, it has been found that the measured S p coefficients are in good agreement with the PT results 
up to p = 7 as long as the variance (<5 2 ) 1//2 is lower than unity 0. This constraint does not seem to 
become tighter for higher orders. 

4.2 From Cumulants to PDF-s 

The shape of the PDF and the values of the cumulants are obviously related. When a limited number 
of cumulants is known it is possible to reconstruct some aspect of the shape of the PDF using the 
so-called Edgeworth expansion (see (2^, 0). But actually as the whole series of the cumulants is 
known for top-hat filtering it is possible to invert the problem (at least numerically) and build the 
PDF from the generating function of the cumulants (see Q for the general method and ||, || for the 
application to PT results). This method is justified by the fact that the S p coefficients converge to 
their asymptotic PT values at roughly the same rate, i.e. for the same values of (<5 2 ) (Fig. 2). 
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Figure 3. Comparison of PT predictions with results of numerical simulations with CDM initial 
conditions for the shape of the one-point PDF of the density contrast. The density PDF has been 

measured at three different radius, R = 5 h _1 Mpc (circles), R = 10 h _1 Mpc (squares), R = 15 
h _1 Mpc (triangles) and compared to the PT predictions calculated from the measured values of the 

variance and the local index (figure taken from ||). 

In Fig. 3, I present a comparison of the PT results obtained in such a way to the numerical 
measurements of PDF-s. Three different smoothing radii have been chosen. For each radius one can 
get the value of the initial index, and the value of the variance. The predicted PDF-s are entirely 
determined by those two numbers. As it can be seen the agreement is extremely good. 

5 Comparison with data 

These successes have given confidence in the validity of the PT results and have boosted investigations 
in the available galaxy catalogues. Using the IRAS galaxy catalogue Bouchet et al. measured the 
3 rd and 4 th moment of the one-point PDF of the local galaxy density. The observed relation between 
the third moment and the second is the following, 

$ RAS ) « 1.5 <^ RAS ) 2 . (21) 

It is a strong indication in favor of Gaussian initial conditions since the exponent is indeed the one 
expected from PT. The coefficient however is lower from what is expected, but the comparison with 
quantitative predictions is complicated for two reasons. The galaxy positions are known only in 
redshift space (that is that their distance is assumed to be proportional to their line-of-sight velocity), 
and IRAS galaxies might be strongly biased with respect to the mass distribution. The first problem 
has been addressed in [^4|] where it is shown with adapted PT calculations that the effects of redshift 
space distortion on S3 are small at large scale. It has been also investigated numerically in |3(J. The 
fact that the measured coefficient, 1.5, is smaller than what is expected (by a factor of about 1.7) is 
then probably due to biases in the galaxy distribution. In case of the IRAS galaxies this is not too 
surprising since those galaxies are known to be under-populated (compared to bright galaxies) in very 
dense areas. 




For these reasons a lot of interest has been devoted recently to the APM angular galaxy catalogue. 
The fact that it has more than 1,300,000 objects makes it the largest galaxy catalogue now available and 
a perfect domain of investigation. Measurements of the S p coefficients have been made by Gaztahaga 
| |19[ | for p < 7. Comparison with PT predictions are however not straightforward because of the 
projections effects with which the relations ( |l8| , |l9| ) for S3 and S4 are not valid anymore. 

Adapted calculations that take into account this new geometry and using the small angle approx- 
imation have been made in f| for any order of cumulants but assuming a power law spectrum. These 
results have been recently extended in for any shape of power spectrum. In fig 4. I present the 
comparison of the measured S3 and 54 coefficients as a function of the smoothing angle (triangles) 
compared to the predicted ones from PT (solid lines). The latter have been calculated with either the 
small angle approximation or with a direct Monte-Carlo integration for which no such approximation 
is required (for S3 only). 
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Figure 4. The s3(9) and sA{9) coefficients in the APM angular catalogue (triangles, see fl9|| ) 
compared to the PT predictions, with the small angle approximation (solid lines), or with a 
numerical integration of the projection effects (squares). 



There is still a significant discrepancy between the PT results and the observations. As the 
theoretical predictions are very robust with respect to uncertainties in the selection function or with 
respect to values of VL and A (see M), it does not seem possible to reconcile the PT predictions with 
the observations, unless galaxies are biased. This bias is however less important than the one expected 
for the IRAS galaxies. These observations are anyway a strong indication in favor of Gaussian initial 
conditions. The implications of models with topological defects like string theories or textures have 
been investigated in particular by Gaztahaga |2(|. He found that the current data, unless the galactic 
distribution is unrelated to the mass distribution^), exclude such theories. 



6 The Velocity Field Statistics 

Although data on the peculiar velocities are more difficult to obtain, it is a very interesting domain 
of investigation since cosmic velocities are in principle directly related to the mass fluctuations. In a 
recent review, Dekel [15] has presented the observational and theoretical status of this rapidly evolving 
field. 

The line-of-sight peculiar velocities are estimated from the Tully-Fisher (or similar) relation be- 
tween the absolute luminosity of the galaxies and their internal velocity dispersion for elliptical, or their 
circular velocity for spirals. From these informations, and taking advantage of the expected absence 
of vorticity, it is possible to build the whole 3D smoothed velocity field. It has been done in particular 



in [11, 17], using the so-called POTENT reconstruction method. A straightforward application of the 



reconstructed velocity field is the use of the velocity-density relationship obtained in the linear regime 
(|6|) that would indeed provide a way to measure Q. This is possible however only if galaxies are not 
biased. Otherwise, assuming that at large scale, ^i a j ax ; es (x) — b 5^(x), the observations constrain a 
combination of Q and b, j3 = Q 0,6 /6. Various results for f3 have been given in the literature (see [|i~5f ). 
A rough compilation of them leads to 0.3 <(3 <1.2. 

As we have no robust models for the large-scale bias of the galaxies, it is quite natural to explore 
the intrinsic properties of the velocity field. In the following, no assumptions are made on the bias, 
galaxies are simply assumed to act as test particles for the large-scale flows. Within this scheme Dekel 



& Rees [16] proposed to use the maximum expanding void to put constraints on Q; Nusser & Dekel 
1 34] tried to reconstruct the initial density field using the Zel'dovich approximation, thus constraining 
on the basis of Gaussian initial conditions. Here I present a more systematic study of the expected 
properties of the local divergence using PT. 

In a similar way than for the density field it is indeed possible to compute the coefficients T p that 
relate the high order moments of the local divergence to the second moment, 

T P = {0V) c /(0 2 ) p - 1 . (22) 

Unlike the S p coefficients these ones are found to be strongly Q dependent (but weakly A dependent), 

T p (U) (x n-°- 6(p - 2 \ (23) 

as is a direct consequence of the time dependence of found in (^). The parameter Tz(£l) was 
calculated in || as a function of 0, 



ft - 1 



f + (n + 3) 



(24) 



and proposed as a possible indicator to measure O. A preliminary investigation using the POTENT 
data gave O > 0.3 with a good confidence level. 
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Figure 5. Comparison of the PT predictions for the shape of the PDF of the velocity divergence 
(solid lines) with results of numerical simulations for n = — 1. In the left panel O = 0.2 and in the 
right = 1. In both cases we have ag ~ 0.4. The dashed lines are the PT predictions when the 

reverse assumption is made on Q. 



But more generally the O dependence of the T p coefficients obviously extends to the shape of the 
PDF of the local divergence. In particular we can show that for a power law spectrum with n pa — 1, 



it is approximately given by, 
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(25) 
(26) 



where cr is the rms of the fluctuations of 0. This distribution exhibits a strong f2 dependence as it can 
be seen on Fig. 5. As expected, when is low the distribution is strongly non-Gaussian, with a sharp 
cut-off in the large divergence tails (it corresponds to rapidly expanding voids). This is this feature 
(discovered in an independent way) that was used in (l6| to constrain Q. 

We have checked these theoretical predictions with numerical experiments (see Fig. 5) and found 
once again that the numerical results are in excellent agreement with them. It opens ways to have 
reliable measure of f2 from velocity data. 



7 Conclusions 

As this rapid overview have shown it, the study of the quasi-linear regime is in rapid development. 
A lot of efforts have been devoted to analytic calculations, and our understanding of the growth 
of structures in the intermediate regime has considerably improved. In particular comparisons with 
numerical simulations have shown that the PT results have a surprisingly large validity domain for 
the density as well as for the velocity fields. So far, most of the comparisons have been done for the 
cumulants at their leading order. However, the recent analytic results obtained for the next-to-leading 
order have open puzzling questions for their interpretation. 

In any case, these PT results provide extremely discriminatory tools to test the gravitational in- 
stability scenarios. Comparisons with data have already provided valuable insights into the properties 
of the galaxy distribution: they have indeed given strong indications in favor of Gaussian initial con- 
ditions and have pointed out precious indications on the existence and nature of biases between the 
large-scale matter distribution and the galaxy distribution. Moreover, assuming Gaussian initial con- 
ditions, it seems possible to get reliable and bias independent constraints on ft from the statistics of 
the large-scale cosmic flows. 
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Resume 

La comprehension de la formation des grandes structures requiert la resolution d'equations non- 
lineaires couplees decrivant revolution des champs de densite et de vitesse cosmologiques. C'est 
un probleme complique qui, ces dix dernieres annees, a ete traite essentiellement avec des simulations 
numeriques a N corps. II y a cependant un regime, le regime dit quasi-lineaire, pour lequel les fluc- 
tuations relatives de densite sont inferieures a l'unite en moyenne. II est alors possible d'utiliser des 
techniques de theorie des perturbations ou les developpements perturbatifs sont faits par rapport aux 
fluctuations initiales. Je presente ici les resultats majeurs qui ont ete obtenus dans ce regime. 



